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ABSTRACT 

Supermassive black holes presumably grow through numerous mergers throughout cosmic time. 
During each merger, supermassive black hole binaries are surrounded by a circumbinary accretion 
disk that imposes a significant (~ 10 4 G for a binary of 10 s M©) magnetic field. The motion 
of the binary through that field will convert the field energy to Poynting flux, with a luminosity 
~ 10 43 ergs _1 (-B/lO 4 G) 2 (Af/lO 8 M0) 2 , some of which may emerge as synchrotron emission at fre- 
quencies near 1 GHz where current and planned wide-field radio surveys will operate. We find that 
the short timescales of many mergers will limit their detectability with most planned blind surveys 
to < 1 per year over the whole sky, independent of the details of the emission process and flux dis- 
tribution. Including an optimistic estimate for the radio flux makes detection even less likely, with 
< 0.1 mergers per year over the whole sky. However, wide-field radio instruments may be able to 
localize systems identified in advance of merger by gravitational waves. Further, radio surveys may be 
able to detect the weaker emission produced by the binary's motion as it is modulated by spin-orbit 
precession and inspiral well in advance of merger. 

Subject headings: black hole physics — cosmology: observations — radio continuum: general — surveys 
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1. INTRODUCTION 

Based on relativistic simulations incor- 
porating force- free elect romag netic fields, 

iPalenzuela. Lehner. fc Lieblingl (j2010f ) suggest that 
mergers of supermassive black holes (SMBHs) which 
occur in the presence of an accretion disk may have 
significant Poynting flux. This Poynting flux may be 
detectable as an electromagnetic (EM) counterpart to 
the gravitational wave (GW) signature of the merger 
(other mechanisms have been proposed as direct and 
indirect electrom agnetic signatures of merger; see 
lSchnittmanll20lil and references therein). These mergers 
will also produce GW signatures, accessible to the Laser 
Interferometer Space Antenna (LISA) (for BH masses M 
in the range ~ [10 3 , 10 7 ]M©) and for exceptionally low 
masses to ground-based GW d etectors (M < 1O 3 M ; 
see, e.g., iReisswig et al.l I2009D . Whether measured 
via GW or EM, the measured merger history will 
strongly constrain our understanding of t he formation 
and evolution of supermassive black hol es (jSesana et al.l 
[27mt ISesana. Volonteri. fc Haardtll2007h . 

Only recently have radio surveys moved beyond inho- 
mogeneous archival data s ets to systematic examinations 
of the variable sky (e .g., iLenc et all 120081 : I Croft et all 
1201 It lOfek et all 1 2 1 lh , and the situation will continue 
to improve. Advances in receivers and digital processing 
make instantaneous fields-of-view of > 10 deg 2 possible 



at GHz frequencies, enabling repeated surveys of wide ar- 
eas of the sky. These technologies are being implemented 
as part of Square Kilometer Array pathfinders unde r con- 
struction ([Johnston et al.ll2007t iBooth et al.ll2009t ). 

While all searches for comp act object merger s 
have so far been negative (e.g., lAbadie et al.l l2010f> . 
the improving performance of both gravi t ationa l (see 
lHarrv fc the LIGO Scientific Collaboration! I2010D and 
electromagnetic surveys increases the discovery poten- 
tial for a wide range of events. In this Letter we consider 
the detectability of the merger event with radio surveys 
centered near frequencies of 1 GHz. We show that the 
flare itself is very unlikely to be detected in the current 
generation of radio surveys, largely independent of the 
amount of EM flux emitted. However, prior to the flare 
there could be other modulation present which may be 
detectable. In what follows, we use a flat ACDM cosmol- 
ogy with S7m = 0.27 and h = 0.72. 

2. ELECTROMAGNETIC COUNTERPARTS OF MERGER 
FLARES 
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IPalenzuela et al.l (|2010t ) simulated the merger of two 
10 s Mq BHs. They found a flare of Poynting flux (with 
I~ 4 x 10 43 erg/s over w 5 hours) that occurred at the 
same time as the GW emission peaked. They also found 
lower-level emission before the flare (L ~ 10 43 ergs _1 ). 
iNeilsen et al.l (|2010l) interpreted the pre-merger secular 
emission as two steady jets powered by the motion of 
^jj each black hole through the background magnetic field, 
with luminosity oc (v / c) 2 B 2 M 2 . For unequal masses, we 
physically expect the luminosity to be provided by the 
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faster, smaller black hole moving through the m agnetic 
field. Using the model of iNeilsen etaU (|2010[ ). if the 
more massive black hole has mass M and the less massive 
has mass qM, we expect a luminosity L oc q 2 M 2 . 

The choice of magnetic field directly affects the elec- 
tromagnetic lu minosities inferred from these simulations. 
Conservatively. IPalenzuela et al.l (|2010l) chose a magnetic 
field that limited their jet luminosity to a small fraction 
L ~ 0.002LEdd of the Eddington luminosity at merger. 1 
We adopt the same assumption: a jet luminosity limited 
to a small fraction eEdd = 0.002 of the Eddington limit 
at the merger event: 

-kflarc = EEdd^Edd (1) 

for q = 1, while for other mass ratios we assume that 
L oc q 2 . 

Given the expected range of magnetic fields (B = 
6 x 10 4 (M/10 8 Mq)" 1 / 2 G, with black hole masses go- 
ing from 10 3 M Q to 1O 1O M0), electrons advected with 
the flow migh t emit synchrotron radia tion near 1 GHz, as 
mentioned bv iPalenzuela et all (|2010l ). Thus, the merger 
flare could be a distinctive radio signature out to cosmo- 
logically significant distances: 



di,Edd =^ \J L/ AttF u 



14.2GpcJg 2 (M/lO 6 M ) 



eradioCeEdd/O-OOS),^ 
(F I/ , min /mJy)(i//GH! 



(corresponding to z f=a 2) where for simplicity we as- 
sume F v oc Ljv. In this expression, rather than model 
the emission mechanism (i.e., spectrum, beaming) in de- 
tail, we assume a fraction e ra dio of this energy is emit- 
ted isotropically in radio frequencies. Given the modest 
Lorentz factor and magnetic field, beaming is not likely 
to be too strong. As much as possible in what follows, 
we attempt to give results that are independent of e ra di - 
Though the emission spectrum is uncertain, the emis- 
sion duration is not: it scales with the total mass of 
the system, as it depends on the orbital t imescale near 
merger. Based on IPalenzuela et all (j2010f ). we estimate 
the merger flare duration by 
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We adopt this estimate for all mass ratios, since the or- 
bital (and hence merger) timescale is determined by the 
more massive BH. 

2.1. Merger Rates 

To assess the visibility of merger flares, we employ 
a merger rate distribution that depends on black hole 
masses and rcdshift. As each comparable-mass merger 
doubles the black hole mass, given the masses and growth 
timescalcs over which they assemble, the supermassive 
black hole merger rate must be both low, less than 

1 The limiting magnetic field required to reach this luminos- 
ity (B ~ 6 X 10 4 G(A//10 8 Mq)- 1 / 2 ) is substantially smaller 
than the magnetic field created by the magneto-rotational insta- 
bility (MRI) at the inner edge of the circumbinary disk, which 
we estimate to be ~ 10 6 G (aM/10 8 M w )~ 7 / 20 l|Pessah et al.ll2006l ; 
Bcgcl man &; Pringlc 20071). Rather than adopt this large circumbi- 
nary field, wc implicitly absorb uncertainties into the ill-determined 
efficiency e Edd . 



10~ 8 Mpc -3 yr _1 , and strongly biased towards low-mass 
mergers: only a few merger events occur per year on our 
past light cone. 

The assembly of SMBHs is reconstructed through 
Monte-Carlo merger simulations, following the hier- 
archical structure formation paradigm. These models 
evolve the BH population starting from BH "seeds," 
through accretion episodes triggered by galaxy mergers, 
and include the dynamical evolution of SMBH-SMBH 
binaries. The SMBH population is consistent with 
observational constraints, e.g., the luminosity function 
of quasars at 1 < z < 6, the M — a relation and the 
BH m ass d ensity at z = (IVolonteri. Haardt. fe Madaul 
20031: IVolonteri. Lodato. fe Nataraianl 120081: 



Volonteri fe Begelmanl 120101) . We adopt two of th e 



fiducial merger distributions used in lArun et all (|2009| ) : 
models LE and SE, where S versus L refers to the 
seed size - lar ge or small - and E refers to "efficient" 
accretion; see ISesana et al.l 1201 lL These models are 
representative of a range of plausible SMBH growth 
scenarios. As with uncertainties in e ra di , we attempt 
to make our conclusions robust to specific merger 
assumptions. 

3. THE VISIBILITY OF MERGER FLARES 

Figure [T] shows the total merger rate as a function 
of BH mass, integrating over redshifts 0-10. Only a 
few mergers per year are expected, even from low-mass 
(< 10 6 Mq) systems. This rate is relevant to untriggercd 
searches by all-sky detectors such as GW observatories 
(LIGO, LISA), which survey the entire sky with roughly 
uniform sensitivity at high duty cycle. For simplicity, in 
what follows we will provide quantitative results primar- 
ily for the LE model; results from the two models are 
comparable for the purposes of this discussion 

For limited-aperture surveys, other factors limit the 
detectable fraction of events (se e the discussion in 
iCordes. Lazio. fc McLaughlinll200"4l for example) . Ignor- 
ing any flux limits, two effects are important. First, sur- 
veys only cover a fraction of the sky CI/A-k, with smaller 
coverage in each pointing. For example, the curvature of 
the Earth restricts telescopes at temperate latitudes to 
Q/4-7T < 80%; individual surveys will cover less. 

Second, surveys return to the same area of the sky 
with a specific cadence T. A telescope which surveys 
a single area continuously (i.e., field-of-view ACl = CI) 
has T ~ 0. More commonly AO <C Cl and the telescope 
spends time doing other tasks. For instance, a survey 
might cover O = 10, 000 deg 2 with 333 pointings of ACl = 
30deg 2 , each lasting 30 s. The survey finishes in < 3hr 
(the smallest possible cadence). If the survey returns 
to each individual pointing 24 hr later, the cadence is 
T = 24 hr. 

With such a survey, the fraction of events that can 
be detected is the fraction that happen to occur when 
observations are ongoing: min[r nare (l + z)/T, 1], assum- 
ing T nar e is much longer than both each pointing and 
any dispersive delay across the bandpass (see § 14. ip and 
simplifying the fla re emission as either on or off (e.g., 
iCenko et all l20ll . In Figure Q] we illustrate how this 
simple cadence cutoff reduces the fraction of low-mass 
merger flares that could be found on the past light cone 
of a survey with cadences T = 1 day, 1 hour, and 10 sec- 
onds. Though many low-mass mergers should occur, the 
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Fig. 1. — Rate of mergers per year on our light cone, for a range 
of total system masses. The thick curves include all mergers for 
two models from[Arun ct al. (20Q3), with solid lines based on the 
"LE" model and dashed lines based on the "SE" model. The thin 
curves reduce this number by max(Tg arc /T, 1) for cadences T = 
1 day (blue), 1 hour (green), and 10 sec (red), assuming the flare 
event duration Tfl are ~ 5h(M/10 8 Mq). 

short durations of their merger flares makes them nearly 
impossible to identify. 

In the LE (SE) model, the total number of mergers 
on our past light cone (summing over all rcdshifts, mass 
ratios, and masses) is about 24 (39) per year. A sur- 
vey with T = 10 s cadence would recover most of the 
events, as even the short low-mass events are sufficiently 
stretched by cosmology that they would be visible for 
AI > 1O 4 M . However, a lh survey is only expected 
to see 3 (2) mergers per year out to z = 10 over the 
whole sky; surveys with finite area will see correspond- 
ingly fewer. Restricting the cadence to 1 d reduces the 
accessible number further, to 0.5 yr -1 (0.2 yr -1 ). 

3.1. Flux Distribution 

So far, we have only counted the number of mergers 
on the light con e of our survey. Using t he predicted lu- 
minosities from iPalenzuela et al.l (|2010[) , we provide in 
Figure [2] the cumulative rate of events greater than a 
flux threshold for a limiting cadence of 1 s (to make the 
numerous low- mass mergers visible). 

Figure [2] shows that merger flares are rare events; the 
flux density corresponding to > lyr -1 over the whole 
sky is only 0.13 mjy (assuming efficiencies of e ra dio = 1 
and CEdd = 0.002) 2 . The brightest events (tens of mjy) 
are extremely rare and come from high-mass mergers, 
generally at moderately high rcdshift. 

In this figure, we have also limited the maximum possi- 
ble duration of a merger flare to 10 6 s (roughly 12 days). 
Longer events are both rare - excluding them changes lit- 
tle - and will be increasingly difficult to localize in time 
and separate from systematic trends. 

4. DISCUSSION & CONCLUSIONS 

Because high-mass mergers are rare (though long- 
lasting) and low-mass mergers produce short and faint 

2 We assume approximately isotropic emission. If the emission is 
tightly beamed, single events will be detectable further away, but 
fewer events will be visible on our lightcone. Since the number and 
timescale of events on our lightcone most limit potential surveys, 
strong beaming will reduce the numbers considered here. 



flares (though common), the rate of potentially de- 
tectable merger flares is small. Even with optimistic 
choices for the efficiencies e ra dio and £Edd, we expect 

< 1 merger per year with the surveys to be conducted 
in the next decade, consistent with zero detections to 
date. Greater flux sensitivity will not increase the de- 
tectable rate substantially, as the finite numbers and 
short timescales limit detectability. Our results depend 
only weakly on the assumed merger rate: while the 
low-mass and low-redshift merger rates are weakly con- 
strained observationally, their merger flares will rarely be 
visible. 

In the radio, ongoing and planned wide-field sur- 
veys have instan t aneous fields-of-view of 1 — 30 deg 2 
(e.g., iCroft et al.l I20I0I : Uohnston et all 120071) at GHz 
frequencies (this increases to several hundreds or even 
1000 deg 2 at a few hundred MHz). Some have rela- 
tively frequent sampling, and cover the same area of 
the sky on timescales from minutes to months, but 
generally only cover a total of < 10 3 deg 2 . Surveys 
that cover a wider area will likely have a cadence of 
at least 1 day. None has the combination of a very 
rapid cadence (ideally < 1 min) and very wide sky 
coverage (> 10 4 deg 2 ) that are likely necessary to detect 
a flare blindly, especially with a required sensitivity of 

< 0.01 mjy. Since the instantaneous fields-of-view and 
cadences of planned optical surveys are typically less 
than or comparable to those of radio surveys and the 
cadence considerations are independent of wavelength, 
optical surveys will be unlikely to discover events like 
these. Only at X-ray and 7-ray energies would planned 
instrumentation be well-suited to the timescales and 
rates of merger flares, although the low fluxes (~ 8 x 

10- 5 (M/10 8 M Q )(^p ho ton/f0kcV)- 1 photon cm" 2 s" 1 
at a redshift of 0.1) might require a next-generation 
mission to be detectable. 

4.1. GW Counterparts 

While Figure [2] suggests that radio flares associated 
with mergers will be difficult to detect, next-generation 
surveys may reach limits more amenable to detections. 
We should consider how merger flares could be identified 
as such and what physics may be learned from them. 

Unlike many proposed counterparts to SMBH mergers, 
this prompt emission mechanism might allow coincident 
detection of electromagnetic and GW signals from the 
same event, even though the circumbinary disk is evacu- 
ated and no accretion onto the compact objects takes 
place. Spatial and temporal coincidence can confirm 
that a radio transient is indeed t he signature of a binar y 
SMBH merger. As reviewed in iBloom fc et all ()2009| ). 
coincident electromagnetic and gravitational signals pro- 
vide an independent cosmological distance ladder, if ac - 
ccssible at cosmological distances (jHofz fc Hughesf2005| ) . 
Additionally, nearby EM counterparts might be local- 
ized to individual host galaxies, allo wing study of galaxy- 
SMBH relations (|Schnittmanll2011l ). 

In contrast with electromagnetic surveys, GW detec- 
tors have roughly uniform all-sky sensitivity at all times 
and a signal that is visible long before merger. Using 
the GW signal as a trigger, electromagnetic observa- 
tions would be freed of the need to survey the whole 
sky continuously; followup observations would be lim- 
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Fig. 2. — Cumulative rate of mergers per year (based on the 
"LE" model) brighter than a given flux density F v in mjy, for 
logarithmic bins of observed duration r = Tfl arc (l + z) in seconds. 
The thin lines are for each bin, while the thick line is for the total 
considering durations of 10° s to 10 6 s. We assumed timescales for 
the flare event that scale as Tfl aro pa 5h(M /10 8 Mq), along with 
efficiency e ra< ji = 1.0 and a frequency v = 1 GHz. 

ited by flux thresholds alone. Here, the large pointing 
uncertainties on current-generation GW facilities (ide- 
ally - 100 deg 2 for a < 10 3 Afo SMB H merger with 
signal-to-noise of 8; e.g.. lFairhurstj|2009D will m ake opti- 
cal follow-up difficult (e.g.. lHaiman et alll2009t) . but arc 
well suit ed to the fields-of-vie w of instruments such as 
ASKAP (j Johnston et al.ll2007| ). Moreover, the GW sig- 
nal may allow identification of an impending merger well 
in advance. For prime LISA-scale sources (10 5 — 1O 7 M ), 
the sky location of an inspiralling binary can be located 
to within 10 deg 2 hours to weeks b e fore the merger even t 
(jMenou. Haiman. fc Kocsisl 120081: iKocsis et all 12007( 1 . 
For the merger trees discussed in this paper, this trans- 
lates to several events per year (slightly less than 1 /3 of 
all £7 &4-detectable even ts) that can be localized this pre- 
cisely ([Arun et al.ll2009 (l . For an optimistic conversion of 
electromagnetic to radio energy, followup pointings will 
identify all LISA events with a flare only if they reach 
a flux sensitivity 0.01 mJy(eEdd/0.002)e ra dio. Less sen- 
sitive followup observations will recover only a fraction 
of events: roughly ~ 17[1 — 0A(logF Vim i n / mjy — 0.5)] 
events per year for F„ im - m <E 0.01 — 3mJy and our fidu- 
cial efficiencies, including all mass ratios. 

For lower-mass mergers (M ~ 10 3 Mq), ground-based 
GW detectors will not identify a sky location before the 
GW merger signal. Nonetheless, if quickly processed, 
their sky localization can still help target EM followup, 
as dispersion delays the EM signal. Plasma dispersion 
will occur in the Milky Way, in the host galaxy of the 
SMBH, and along the line-of-sight in the intergalactic 
medium. For cosmological sources the total dispersion 
measure (DM, the integral of th e electron column 
density) may reach > 1000 pc cm - 3 (jlnouel 12004( 1. This 
implies a time delay At = 4.15^q 2 Iz DM ms, requiring 
rapid localization and repointing. Even for short flares 
with modest dispersions, such a delay would be hard 
to detect, but it is possible if the radio cadence is 
sufficiently short or the observing frequency low. It 
may also be possible to detect dispersion within the 
radio data itself, by measuring the relative delays of 



different frequencies across a bandpass of width Av 
(St = 8^Qg z A^GH Z DMms). This is more difficult, 
since across a finite bandpass the rel ative delay is 
even smaller, but is routinely done (e.g.. lLorimer et al.l 
12001 . In fact, for very low mass events M < 
lO 4 M (DM/lOOOpccm- 3 )(^/lGHz)- 3 (A^/3OOMHz) 
dispersive smearing will exceed Tfl arc , but this will 
not greatly change Figure [TJ However, matching this 
"internal" delay with that relative to GW observations 
could prove a powerful confirmation of the nature of the 
event. 

Second-generation ground-based GW detectors will be 
sensitive to the lowest-mass mergers (M ~ 200 — 10 3 Mq) 
out to a strongly mass- and orientation- dependent 
threshold z ~ 0.1 — 2. The associated EM flares will 
be short (dispersion-limited) and faint. With the most 
optimistic efficiencies eEdd) Cradio, radio surveys would 
have comparable reach to GW surveys at F v ~ 0.1 mjy 
(Eqn. [2]); with less efficient conversion or followup, fewer 
coincident events can be found. Unfortunately, unlike 
SMBH mergers, observations do not directly constrain 
such merging binaries. The low-mass mergers to which 
ground-based detectors are sensitive simply may not oc- 
cur. Even if they do, both GW and radio observations 
are sensitive to a minute fraction of t he universe (not tru e 
for third-generation GW detectors; ISesana et al.l l200"9( l . 
That said, if radio surveys can distinguish short (< 1 s) 
flares in targeted observations, ground-based GW detec- 
tors working in concert with radio telescopes can rule 
out extremely optimistic (> 10~ 8 Mpc" 3 yr _1 ) low-mass 
SMBH merger rates and efficiencies. 

4.2. Non-Merger Events 

While the rate of potentially detectable mergers is 
small, other EM emission associated with binary SMBH 
inspiral could be detectable. EM emission from SMBHs 
is well known across a range of wavelengths; radio emis- 
sion from active galactic nuclei is common. Wc differ- 
entiate between generic AGN emission and that associ- 
ated with an orbiting pair of SMBH through the time 
domain. AGN do vary intrinsically but mostly ape- 
riodically; detecting such periodic behavior in a radio 
light curve would be a strong indicat ion of an inspi- 
ralling SMBH pair (e.g.. lKomossall2006( l: we defer addi- 
tional methods for confirmation to a forthcoming paper 
(O'Shaughnessy et al. 2011, in prep.). A number of bi- 
nary AGN are known or suspected (e.g., iKomossal I2QQ& 
Rodriguez et alJl2006t [Smith et al J 120101: iBurke-Spolaorl 
20101) . Most of these have evidence from a resolved pair 
of bright spots or a double set of emission lines, but they 
all pr obe systems far from actual merger (jBurke-Spolaorl 
|2010D . We consider what might happen as the systems 
approach merger. 

Variability will happen over a range of timescales. 
First, even before the merger the E M flux of the sys- 
tem is expected to increase as (v/c) 2 (lMcWilliamsJ(2Tjl0l: 
iNeilsen et al.ll2010l ). where v ~ (imerge — t)^ 1 / 8 traces the 
increasing orbital speed during inspiral, going to a max- 
imum of u max ~ c/V6 « 0.4c at the innermost stable 
circular orbit, and with a singularity at merger (t merge ). 
The flux increase will be secular and may be detectable, 
but the timescales over which it changes appreciably are 
likely either too long (during the lengthy inspiral) or too 
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short (right before merger), and will be difficult to iden- 
tify uniquely. 

A promising ca ndidate is v ariability induced by pre- 
cession (also see iKatzl [19971 for a related discussion). 
If there is a Poynting flux associated with a jet, as in 
iPalenzuela et al.l ()2010t ). the axis of this jet could precess 
if the BH spins arc not aligned with the orbital angular 
momentum. This would presumably cause the EM sig- 
nal to vary o n that time scale (although it could be more 
complicated; lKatzlll997T) . The precession timescale is ex- 
pected to be t p ~ M(y/c)~ b (jApostolatos et al.lll994l ). 
The scale of the variations is not known (it depends 
on the anisotropy of the emission), but could easily be 
> 50%. 

While a full treatment is beyond the scope of this pa- 
per, we are drawn to consider the detectability of pro- 
cessing jets in binary SMBHs for two reasons. First, 
the time spent at a moderate velocity v/c ~ 0.1 com- 
pared to the duration of the merger itself is large, scal- 
ing as (v/wmax) -8 - A much larger number of systems 



exist in this state compared to those merging; their 
timescales are more amenable to detection. Second, a 
system will undergo many precession cycles, so periodic 
modulation may be detectable (along with other changes, 
such as secular increase or orbital modulation) ; we expect 
Ap ~ (v/v max )~ 3 periods to be visible in a roughly log- 
arithmic velocity range. Precession has likely been seen 
in galactic BH binaries (|Katz|[T997t) . and does have ob- 
servational consequences for the jet emission. In a future 
paper (O'Shaughnessy et al., in prep), we will discuss the 
detectability of binary SMBH jet precession in detail. 
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